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SUMMARY 


Sonic-boom wave shapes and caustic lines generated by an airplane per- 
forming a general maneuver are studied. The equations are programmed for 
graphical output as a perspective view of the wave shape. This quasi-three- 
dimensional presentation provides a qualitative insight into the effects of 
the maneuver on the wave shape and the caustic locations. 

For the special case of planar maneuvers , the principal curvatures of 
the wave front are derived. These curvatures are needed to calculate the 
sound field in the vicinity of a caustic. 

The results of the analysis are applicable not only to sonic-boom studies 
but also to the calculation of noise generated by a supersonic rotor or pro- 
peller blade tip. 


INTRODUCTION 

When a supersonic aircraft accelerates or turns, the wave front generated 
by the flight trajectory develops folds or caustic lines. Special interest is 
attached to the regions near caustics because the sonic boom overpressures are 
enhanced in such regions. The sonic-boom literature contains a number of anal- 
yses which treat the local conditions near a caustic. (See, for example, 
refs. 1 to 3.) Heretofore, these analyses have assumed expressions for the 
local wave-front shape and curvature which are required to perform the calcula- 
tion in the neighborhood of the caustic. The problem that remains is to relate 
the wave-front shape and curvature to the actual maneuver that is being per- 
formed by the aircraft generating the wave. The solution to this problem is 
the primary subject of this paper. 

Because of the complexity of the analysis, the expressions derived from 
the wave-front curvatures are limited to planar maneuvers. This means simply 
that the torsion of the flight trajectory is zero. The maneuver plane may be 
horizontal, vertical, or skewed. 

On the other hand, the equation for the rays and the wave-front shape is 
not limited to planar maneuvers. Some types of sonic-boom problems, especially 
those that are strongly dependent on refractive effects, are more advantageously 
formulated in terms of the rays (refs. 4 and 5) . However, when the flight pro- 
gram involves turns and acceleration, the wave shape is no longer obvious, and 
it is desirable to know the location of the entire wave front at any given time. 
Other advantages of the wave-front formulation are discussed in reference 1 . 
However, from an engineering standpoint, probably the greatest advantage of the 
wave-front formulation is the intuitive understanding of the effects of maneu- 
vers that it can provide. Toward this end, the wave-front equation has been 
programmed with graphical output in the form of a perspective view of the wave. 
This type of quasi-three-dimensional picture provides a rapid qualitative 



insight into the relationship of the wave-front shape and caustic locations 
to the aircraft maneuvers. 


SYMBOLS 

In numeric examples, distances are nondimensionalized in terms of the 
distance that sound travels in 1 second. 

a speed of sound 

b = T x n 

A,B,C,D labels for groups of expressions defined by equations (A4a) to (A4d), 
respectively 

E,F,G coefficients of first fundamental form for a surface, equations (Alb) 
to (Aid), respectively 

e,f,g coefficients of second fundamental form for a surface, equations (A2a) 
to (A2c) , respectively 

i,j,k orthonormal base vectors, fixed in space 

M Mach number 

n unit vector perpendicular to flight direction, in plane of maneuver 

P curvature parameter (see fig. 8) 

R = a (t - T) 

r generic vector of point on shock wave 

r^- position vector of aircraft 

s distance along flight path 

T unit base vector in direction of local tangent to flight path 

t time, sec 

V local flight speed 

3 = \/m 2 - 1 

A denotes an increment in a quantity 

0 angle variable for curved flight path 

K total surface curvature 
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curvature 


K 

P 1 /k radius of curvature of flight path 

PlfP2 principal curvatures of wave front 

T time at which aircraft is at position ? 

$ angle variable in plane of characteristic, measured from horizontal 

Subscripts: 

t flight path or trajectory 

T derivative with respect to T 

<t> derivative with respect to 4> 

1 first 

2 second 

A prime denotes a derivative with respect to time. 


ANALYSIS 

Wave-Front Shape 

It is convenient to use twosystems of unit base vectors. One is the 
usual system of unit vectors (i,j,k) fixed in space, with k taken normal to 
the plane of the_maneuver for planar trajectories. The other system of unit 
base vectors (T,n,b) moves with the aircraft. (See fig. 1.) The first vector, 
T, is in the direction of flight and is defined by 


dE t 

T = — (la) 

ds 


or 


df t 

= VT (lb) 

dT 
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The second vector, n, is defined by 


dT 

— = K t R (2a) 

ds 


or 


dT 

— = VK+n (2b) 

dT 


where K t is the curvature of the flight path. (See ref. 6, p. 18.) Since T 
is a unit vector, n is normal to T in the plane of the trajectory. For planar 
maneuvers, the third vector 5 is parallel to k but may be oppositely directed. 

The flight program may be defined in several different forms. It may be 
given as a vector position of time r t (x), as a trajectory r t (s) and velocity 
function of arc length V(s), or as a trajectory r t (s) and velocity function 
of time V(T). in any case, the time and arc length parameters can be related 
by 


T = 



ds 

v 


or 


s 



V dT 


At time t the pressure disturbance created by the moving aircraft at 
time T has spread a distance R = a(t — T). If r(t,T) is the generic point 
of this disturbance, and f^.(T) is the source location at time t, then the 

equation for the disturbance wave front is 

(r - r t ) • (r - r t ) = a 2 (t - T) 2 = R 2 (3) 
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The shock surface is the envelope of this one-parameter (x) family of 
spheres. The characteristic lines that comprise this envelope are found by 
differentiating equation (3) with respect to T and solving the resulting 
equation simultaneously with equation (3) . Thus, 


(r - r t ) • VT = aR 


(4a) 


or 


(r - r t ) • T = R/M 


(4b) 


Thus the component in the flight direction of a ray to the shock envelope 
is R/M, which means that the characteristic line is a circle with center at 
r^ + (R/M)T and with radius Br/M. In the plane of this circle let 4 1 denote 
the angle with the horizontal. (See fig. 2.) Then the vector equation for the 
shock surface at time t can be written explicitly in terms of the moving 
trihedral base vectors and the two parameters T and $ as follows: 


a (t - T) aB(t) (t - T) cos 4> aB(x) (t - X) sin 4> _ 

r = T(x) + — : h(x) + —7: b(X) + r t (5) 


M (X) 


M (X) 


M(X) 


If, on the other hand, both x and <t> are held constant and t is 
varied, then equation (5) describes the trajectory of an element of the wave 
front, that is, a ray. 

The equation for a caustic line in the wave front at time t is found by 
differentiating equation (4) with respect to x and solving the resulting 
equation simultaneously with equations (3) and (4). This is the mathematical 
equivalent of the geometric condition that two consecutive characteristic cir- 
cles be tangent to each other. The derivative of equation (4) with respect to 
X is 


• VT + (r - r t ) • V 2 — + V'T) = -a^ 
\ ds 1 


dT 


= 


Substituting from equations (1) and (2) yields 


(r - r t ) • (V 2 K t n + V'T) = 3 2 a 2 


( 6 ) 
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Combined with equation (4) , this equation becomes 


/B 2 v'r\ 

(r - r t ) • n = P t — - — (7a) 

\M 2 V 2 M/ 

Equations (4) and (7a) , respectively, give the components in the T and 
h directions of the ray to a point on the caustic. The component in the b 
direction follows from equation (1) as follows: 


/ B 2 r 2 

.ZB 2 v'r\ 


(r r t ) * b ± 1 

- p t 2 , - „ 

(7b) 

\j M 2 

\M 2 V 2 M / 



Equations (4) and (7) were given in reference 7 in a slightly different 
form. In order to calculate the wave shape from equation_(5) or the caustic 
line from equations (4) and (7 ) , the_moving vectors (T,n,b) must be expressed 
in terms of the fixed base vectors (i,j,k) by using equations (1) and (2). 

In the case of a straight accelerating flight, «{. is 0 and the coeffi- 
cient of n in equation (6) becomes zero. For this axisymmetric case, the 
caustic line at time t is the circle determined by the value of x that 
satisfies the equation 


3 2 (t)M(t) 

t + = t 

M'(x) 


which is obtained from equations (4) and (6) . 


Wave-Shape Examples 

Even relatively simple maneuvers present some interesting examples of 
wave shapes. Figure 3 shows in perspective the caustic line and a number of 
the characteristic circles comprising the wave front that results from a simple 
turn at constant speed. Figure 4 shows how the wave folds on itself to form a 
caustic when the flight trajectory changes from a straight flight to a turn. 
Only the half of the wave above the flight plane is shown in order to avoid 
the confusion of too many overlapping lines. Figure 5 shows a different per- 
spective of a similar maneuver , but with the airplane decelerating at a con- 
stant rate from an initial Mach number of 2.4. The principal effect of the 
deceleration is to move the caustic line farther from the flight trajectory 
than it would be for constant flight speed. Figure 6 shows a similar maneuver 
but with the airplane accelerating from an initial Mach number of 1.1. Two 
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caustic lines appear in this case. The first occurs near the beginning of the 
acceleration region and is due to the acceleration, while the other is associ- 
ated with the turning maneuver. 

The nature of the acceleration caustic is better seen in a detailed plot 
of the wave cross section (fig. 7) . At seme distance from the flight path a 
single caustic line forms initially, but two caustic lines appear as the wave 
front crosses itself. These lines appear in the cross-sectional plot as cusp- 
like points of reversal of direction. Note, however, that they do not repre- 
sent cusp or arete points in a caustic line. Some experimental data on sonic- 
boom ground patterns and signature shapes are presented in reference 8. 


Determination of Wave-Front Curvature 

In order to calculate the flow field in the vicinity of a caustic, it is 
necessary to know the curvature of the wave front approaching the caustic line 
(refs. 1 to 3) . In addition, knowledge of the curvature distribution over the 
wave front provides an insight into the relative expansion or compression of 
the wave. 

The shock wave envelope is, except for straight, constant-speed flight, a 
surface of double curvature. At a fixed time t, the equation of the surface 
(eq. (5)) is a function of the two variables T and <)>. The theory for calcu- 
lating the curvatures with the equation in this form is described in chapter 2 
of reference 7. The procedure is to calculate the total radius of curvature, 
which is the reciprocal of the total curvature. The details of the calculation, 
which are somewhat tedious, are given in the appendix. The result is as 
follows: 


1 

K 


Pi P 2 = -a ( t - T) 


8a 


M' 

— + MaK t cos <j> 
8M 


- a(t - T) 


( 8 ) 


Some insight into the nature of the principal curvatures can be gained by 
examining the example of straight, accelerating flight. Since for this case 
K t = 0, then D, f, and F are all zero (see appendix) . Consequently, the 
curvature directions are orthogonal and the principal curvatures can be computed 
by the following equations (see ref. 6, p. 81): 


e M' A M' -8 b M' 

E B 2 A 2 + B 2 B 2 M 2 B 2 BM 2 B 
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Pi 


3m 2 b 

M' 


(9a) 


Therefore, 


Substituting for B gives 


3 2 Ma 

Pi = a(t - t) (9b) 

M' 


and 


G MC 

P2 = “ = 7- = a (t * T > (10) 

g 3 


Thus, tO| is the curvature of the meridians along the shock envelope. 
This curvature results from the acceleration maneuver. On the other hand, 

<2 is the curvature associated with the normal acoustic spreading of the wave 
front. 


Equation (1 0) could therefore have been obtained by the following argument. 
Since, according to equation (4) , the wave surface consists of circles of radius 

3a (t - T) 

, the radius of curvature normal to the surface is this radius divided 

M 

by the cosine of the angle between the normal to the surface and the radius vec- 
tor of the circle. This result is in accordance with Meusnier's theorem (ref. 6, 
p. 76) . Since this cosine is 3/M, the normal radius of curvature is simply 
a(t - T) , which is the same as equation (10). 

This latter argument is applicable locally not only for the case of 
straight flight but also for a general maneuver. At any point on the shock 
surface one of the principal radii of curvature is P 2 = a (t - x) . Consequently, 
according to equation (8), the other radius of curvature is 


3 2 Ma 

Pi = a(t - x) (11) 

M* + 3M 2 a< t cos <|> 


The curvature k-| results from the maneuver. The curvature direction asso- 
ciated with K-| is locally normal to the characteristic circle. Consequently, 
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near a caustic, where p-j =0, the curvature direction is nearly normal to the 
caustic, since the caustic is locally tangent to a characteristic line. 

Thus, in the vicinity of a caustic line whose curvature is small compared 
with K*| , the wave is approximately two-dimensional in character with its 
radius of curvature given by equation (11). 

The essential expression in the curvature calculation is the last term in 
equation (11). This term (denoted by P) represents the distance of a point 
on the caustic line from the trajectory point at which it was emitted. It is 
readily seen that a positive acceleration decreases this distance and the 
radius of curvature at a given distance along the ray, whereas a deceleration 
has the opposite effect. These results are shown in the plots of figures 8(a) 
and 8(b). The Mach numbers (1.2 and 1.6) for which these examples were calcu- 
lated are the instantaneous Mach numbers at the time the caustic ray was emitted. 
Fran the parameter P plotted in figure 8, the curvature P] can be calculated 
as R - P, where R and P are the first and last terms on the right side of 
equation (11). The total curvature is then 1 /Rp-| . Equation (11) could be 
obtained in a brief but intuitive manner as follows: 

Substitute 


U - r t ) 


T = 


a(t - T) 
M 


and 


(f - f t ) 


n 


3a (t - x) cos 4> 
M 


into equation (6) . This yields 


a(t - x) (V 2 K t B cos 4> + V') = B 2 Ma 


(12a) 


or 


a ( t - T) 


B 2 Ma 

= 0 

M 1 + 3M 2 aK t cos 4> 


(12b) 


This is the condition that exists at a caustic point. It can also be obtained 
by setting p-| = 0 in equation (11). Therefore, by using dimensional 
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considerations and the fact that the rays are straight, one can argue that 
equation (11) follows from equation (12) without making the extensive curva- 
ture calculation. 


Ray-Tube Area 

An incidental result of the preceding calculation for the curvature is a 
formula for the ray-tube area. As given on page 63 of reference 6, an elemental 
area of the wave front is determined by the square root of the discriminant of 
the first fundamental form: 

Aa = \J EG - F 2 dT d<J) 


Thus, the variation of the cross-sectional ar ea of a rav tube is deter- 
mined by allowing t to vary in the expression \/eG - F 2 , while holding x 
and 4> constant. Equations (A4a) , (A4c), and (A8) yield R = a (t - X): 


S 2 a 

Aa = — r 

M 


^ I M' M< t 00:3 9 \ 


,B 2 Ma 


B 


dx d<J> 


( 13 ) 


Equation (13) is the same result that was obtained by Rao (refs. 9 and 10) 
by means of a somewhat more involved development. It may also be noted that 
equation (5), for the wave shape, and equation (13), for the ray-tube area, are 
useful not only for sonic-boom calculations associated with maneuvering super- 
sonic airplanes but also for computing the noise field of supersonic propeller 
or rotor tips. 


Higher -Order Caustic Locations 

A caustic point represents a point of tangency of two successive charac- 
teristic circles of the wave front. Under certain conditions, caustic points 
may also exist at which a higher-order contact occurs among successive charac- 
teristics. Additional compression, beyond that normally experienced at a 
caustic point, would be realized at these points. 

The location of such a point is found by differentiating equation (6) 
with respect to x and solving the resulting equation simultaneously with 
equations (3), (4), and (6). The derivative of equation (6) is 

-VT • (V 2 K t n + V'T) + (W'K t h + V"T - V 3 K fc 2 T + 2W'< t n 

+ V 2 K^.n) • (r - r t ) = 2W' 


1 0 



or, with terms collected. 


(r - r t ) • p3W'K t + V 2 Kt)n + (v" - V 3 K t 2 )TJ = 3W' (14) 

For some maneuvers, equation (14) does not have a solution. One example 
is a straight linear acceleration, for which ic t = 0 and V" = 0. Another 
example is a constant speed, constant radius turn, for which V' = 0 and 
K\. = 0 . 

It should also be noted that, even if it is possible to solve equation (14), 
the caustic line does not necessarily have a cusped shape. A simple example is 
straight nonlinear acceleration (V" # 0) . If a value of t - x exists for 
which both equations (6) and (1 4) are satisifed, then there is a higher-order 
folding of the wave along the corresponding characteristic circle, but no cusp 
is formed. 


CONCLUDING REMARKS 

The relationship of sonic-bocm wave shapes and caustic lines to the maneu- 
ver performed by the generating aircraft has been studied. The analysis treated 
general turning and accelerating maneuvers in an isothermal atmosphere. The 
equations were programmed for graphical output in such a way that the wave 
shape could be observed in a quasi-three-dimensional perspective view in order 
to provide a qualitative insight into the effects of the maneuver on the wave 
shape. 

For the special case of planar maneuvers, the principal curvatures of the 
wave front were calculated. These curvatures are of interest in the calcula- 
tion of the sound field in the immediate vicinity of a caustic line. 

Although the analysis is based on equations originally derived for a 
maneuvering supersonic airplane, the results are also applicable to the propa- 
gation of sound from a supersonic rotor or propeller tip. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 7, 1979 
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APPENDIX 


METHOD OF CALCULATION OF TOTAL RADIUS OF CURVATURE 

This appendix provides the details of the calculation of the wave-front 
curvature K, which is the product of the two principal curvatures, k-j and 
i< 2 . The basic formula, which is given on page 83 of reference 6, is 

eg - f 2 

K = (Ala) 

EG - F 2 


In this expression, the numerator is the discriminant of the second funda- 
mental form of the equation of the surface, and the denominator is the discrim- 
inant of the first fundamental form. The individual quantities are the coeffi- 
cients which are defined by 


E = r T • r x 

(Alb) 

F = r x • r<j, 

(Ale) 

G = r<j, ♦ r<j> 

(Aid) 


Also required are the coefficients of the second fundamental form, defined by 
e = (r xx 


r T x r$) j^EG - F 2 

(A2a) 

r x x r<j))/^EG - F 2 

(A2b) 

r x x r cj ,)/ \|eG - F 2 

(A2c) 


The derivatives of r are obtained by differentiating equation (5) 


r T = 


B 2 a M'a(t - T) 


M 


- BKt-a 2 (t - T ) cos $ 




T + 


a 2 K t (t - T) + 


M'a(t - T) cos cj) 

6m 2 


3a cos 
M 


n + 

M'a(t - T) 

Ba" 

Bm 2 

M 



sin 4> b 


(A3) 


1 2 



APPENDIX 


In order to simplify the calculations , the mathematical expressions may be 
grouped as follows: 


B 2 a M'a(t - T) - 

A = B<*.a^(t - x) cos (j) (A4a) 

M M 2 


M'a(t - T) 6a 

B = 

6m 2 « 


(A4b) 


C = 


Ba(t - x) 
M 


(A4c) 


D s K t a 2 (t - T) 


(A4d) 


Thus, 


A = -6(B + D cos <J>) 


(A4e) 


In terms of these groupings, equation (A3) becomes 


r T = AT + (D + B cos <|>)n + B sin <J> b 


(A5) 


The derivative of equation (5) with respect to <j) is 


r<j) = -C sin 4> n + C cos <J> b 


(A6) 


The coefficients of the first fundamental form can now be obtained from equa- 
tions (AT ) , (A5) , and (A6) : 


E = A 2 + D 2 + 2BD cos <J> + B 2 


(A7a) 


F = -CD sin <(> 


(A7b) 



APPENDIX 


G = C 2 (A7c) 

The discriminant of the first fundamental form is found by substituting equa- 
tions (A7) into the expression EG - F 2 and simplifying by using equations (A4) . 
This yields 


EG - F 2 


m 2 c 2 a 2 

B 2 


(A8) 


The second derivatives are required in the calculation of the coefficients of 
the second fundamental form: 


r XT = [A' - VK t (B cos <j> + D)]t + (VK t A + B' cos 4> + D')n + B' sin 4> b (A9a) 


r T( j) = BD sin <j) T - B sin (j> n + B cos <J> n (A9b) 

= -C(cos <}) n + sin 4> b) (A9c) 

From equations (A5) and (A6) , 

x t,j) = C(B + D cos 4>) T - AC cos <|> n - AC sin 4> B 
or, from equation (A4e) , 


r T x r 4> = 



+ cos 4> n + sin 



(A10) 


The first coefficient, e, is obtained from equations (A2a) , (A9a) , (A8) , 
and (Al 0) : 


M 


— e 

B 



- VK t (B cos (() + D)J 


cos 4>(VK t A + B' cos 4> + D') - B' sin 2 4> 


(All) 
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From equation (A4e) , 


MM' 

A' = -8(B' + D' cos (J>) - (B + D COS 4*) 


MM' 

= -3 (B ' + D' cos <j>) + A 

B 2 


Consequently, the first term and the last three terms on the right side of 
equation (All) combine to yield MM'A/3 3 . Thus, 


M'A w t r 

e = + (B - BA) cos cj> + d 

82 M 


However , 


B cos 4> + D = (B + D cos <l>) cos <j> + D sin^ <t> 


A cos <f> , 

+ D sin^ <}> 

8 


Therefore, 


e = 


/m' w t VK fc 8 cos <})\ 

- — + - — cos (j) + A + aKj.D sin 2 $ 

U2 8m M 


or 


m 1 M 2 aK t cos <j>\ 

e = -| — + A + aKf-D sin 2 4> 

[ 6 2 B 1 


(Al 2a) 
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The second coefficient, f, is obtained from equations (A9b) and Al 0) 


/ ^ \ 

f = -ACD sin 4> 

\MCA/ 


f = - 


BD sin <{) 
M 


(Al 2b) 


Similarly, from equations (A9c) and (Al 0) , 


g = AC 2 [— ] 

\ MCA/ 


Be 

g = — 

M 


(Al 2c) 


The discriminant of the second fundamental form can be obtained from 
equations (Al 2) : 


/m* \ BaK t CD sin 2 4> g2 D 2 sin 2 ^ 

eg - f 2 = - — + MaK t cos <J>] CA + — 


,Bm 


M 




Since CaK t = Bd/M, the last two terms cancel: 


o /M’ \ 

eg - f z = -I — + MaK t cos 4>j CA 


(Al 3) 


The total curvature can then be calculated by using equations (A8) and (Al 3) 


K = k-\K 2 = 


eg 


- f 2 


EG - F 2 M 2 AC\^ m 


B 2 /M' \ 

— + MaK t cos cp 
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APPENDIX 


It is more convenient to work with the reciprocal, which after substituting 
from equation (11), is 


- = PlP2 = -a(t - t) 
K 


8a 


M' 


- a(t - t) 


— + MaK*- cos <)) 
8 m fc 


(AT 4) 
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Characteristic 
circle of wave front 


Figure 2.- Ray and wave-front geometry. 



Figure 3.- Perspective view of wave front with caustic resulting from a 
constant-curvature turn at M = 1.16. 
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Normal view of flight path 
(Turn radius = 22.9 km) 

Figure 4,- Portion of wave front above flight plane for a straight 
flight path succeeded by a constant-radius turn at M = 1.3. 
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Normal view of flight path 
(Turn radius = 22.9 km) 


Figure 5.- Portion of wave front above flight plane for a straight flight path 
succeeded by a constant-radius turn with aircraft decelerating from M = 2.4 
to M = 1 ,52 at dM/dT = -0.004/sec, 
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Flight path 



Normal view of flight path 
(Turn radius = 22.9 km) 

Figure 6.- Portion of wave front above flight plane for a straight flight path 
succeeded by a constant radius with aircraft accelerating from M = 1 .1 to 
M = 3.1 at dM/dT = 0.016/sec. 
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Figure 8*- Wave-front curvature parameter P = as a function 

M' + 3M 2 aK t cos <|> 

of trajectory radius of curvature for various values of acceleration. 
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